A REMARK ON THE EQUIVALENCE OF ISOKINETIC AND ISOENERGETIC 
THERMOSTATS IN THE THERMODYNAMIC LIMIT. 



by David Ruelle*. 



Abstract. The Gaussian isokinetic and isoenergetic thermostats 
of Hoover and Evans are formally equivalent as remarked by 
Gallavotti, Rondoni and Cohen. But outside of equilibrium the 
fluctuations are uncontrolled and might break the equivalence. 
q ! We show that equivalence is ensured if we consider an inhnite 

CN ' system assumed to be ergodic under space translations. 
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1. Introduction. 



In the study of nonequilibrium statistical mechanics, if nonhamiltonian forces are 
used to achieve nonequilibrium, a thermostat is needed to cool the system. The Gaussian 
thermostats introduced by W. Hoover and D. Evans have the great interest of respecting 
the deterministic character of the equations of motion (see for instance Evans and Morriss 
[3]). Starting with an evolution equation x = F(x) in phase space, a Gaussian thermostat 
constrains the evolution to a prescribed hypersurface E by projecting F(x), for x G E, to 
the tangent plane to E at x. In the present note we follow Cohen-Rondoni, and Gallavotti 
comparing an isokinetic and an isoenergetic thermostat, and showing that they give the 
same result in the limit of a large system (thermodynamic limit). 

In equilibrium statistical mechanics one can show rigorously that fixing the kinetic 
energy is equivalent to fixing the total energy, asymptotically for large systems (see [6]). It 
is therefore natural to hope that something similar is true for nonequilibrium, as advocated 
by Gallavotti (many references, see [4], [5]) and by Cohen and Rondoni [2]. However, the 
entropy considerations which are available in equilibrium statistical mechanics fail utterly 
outside of equilibrium, i.e., fluctuations of energy at fixed kinetic energy are uncontrolled, 
and the situation appears rather hopeless. We shall show however that the argument of 
Cohen and Rondoni can be modified to apply, at least formally, to the dynamics of actually 
infinite systems. (In a different context - at equilibrium - Sinai [7] has also shown the 
interest of considering the dynamics of infinite systems). Our approach will remain formal 
at the level of infinite system evolution equations: technical problems arise there, which 
do not seem directly related to the problem at hand, and are better discussed separately. 

We shall consider a system of particles in d dimensions which is infinitely extended 
in v dimensions, with 1 < v < d, and we shall discuss states of infinitely many particles 
which are invariant under translations in R^. The assumption that the infinite systems 
dynamics is well defined, and R^-ergodicity, will be sufficient to establish the equivalence 
of isokinetic and isoenergetic nonequilibrium steady states 

2. IK an IE dynamics. 

We recall now the definition of the Gaussian isokinetic (IK) thermostat. We take for 
our configuration space M a compact subset of R u x T" where T" is the w-torus, and 
momentum space is identified with Tt u+V . We assume that a force field on M is given, 
written as — gradF + £, where V : M — > R is a potential, and £ is a nongradient vector 
field*. Consider now the equations of motion 

p = -d q V + i-ap 
q = p/m 

completed by elastic reflection at the boundary of M. Without the term £ — ap this time 
evolution would be Hamiltonian. The term £ maintains the system outside of equilib- 
rium. The term — ap is the thermostat. We obtain the Gausssian isokinetic thermostat by 

* Note that a change in V can be compensated by a corresponding change in £: the 
splitting of the force into two terms is arbitrary for the IK time evolution. 
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choosing a such that the kinetic energy is constant: 



i.e., 

a = (-d q V + 0-p/p 2 (2) 

Note that if £ is locally a gradient (corresponding to a multivalued potential function on 
M), the Dettmann-Morriss pairing theorem asserts that (except for one value = 0) the 
spectrum of Lyapunov exponents of an ergodic measure is symmetric with respect to some 
constant c which is in general nonzero. (We shall however not make use of this result). 

We consider now the Gaussian isoenergetic (IE) thermostat associated again with the 
force — gradF + £, but where we want to maintain fixed the energy function 

H = p 2 /2m + V(q) (3) 

The equations of motion are again of the form (1) and using (3) the isoenergetic condition 
is 

= H=P-- (-d q V + £ - ap) + d q V ■ 

m m 

i.e., 

a = £ • p/p 2 (4) 

With the Gaussian isoenergetic (IE) thermostat the time evolution is thus defined by (1), 
(4). 

We consider now the IK and the IE time evolution in the infinite system limit. We 
want to study the time evolution of a state p ergodic under R^-space translations. We 
shall ignore existence and uniqueness problems for these evolution equations, and our 
discussion will thus remain formal in this respect. (In fact, the one-dimensional situation 
may be relatively accessible to rigorous study, but the n-dimensional case with n > 2 
appears much more difficult). 

Physically we may think of a system of particles in a region D invariant under Tt u , 
where 1 < v < dimD but possibly v < dimD. For example we may consider a shear 
flow between two moving plates, but we do not take the limit where these two plates are 
infinitely far apart, as this would introduce unwanted hydrodynamic instabilities. Another 
example would be a system of particles in [0, L] x Tt u . In the x-direction we put an electric 
field and we assume a suitable boundary condition (see [1]). 

The expressions p 2 = p-p and £-p diverge for an infinite system, but behave additively 
with respect to volume, and we can (under mild conditions on p) define the average per 
unit volume with respect to p, noted {p 2 ) p or (£ • p) p . Since p is ergodic, it is carried by 
points (in infinite phase space) for which the large volume average of p 2 or £ • p is well 
defined and constant, equal to (p 2 ) p or (£ ■ p) p . The expressions V, d q V -p behave almost 
additively with respect to volume and, again under mild conditions, we can define the large 
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volume averages (V) p , (d q V ■ p) p . Again, p is carried by points (in infinite phase space) 
for which the large volume average of V or d q V ■ p is well defined and constant, equal to 
(V) p or{d q V-p) p . 

In our formal treatment of the infinite system IK or IE evolution we consider the 
time evolution of an infinite phase space point, generic with respect to the space ergodic 
measure p, replacing the expressions (2), (4) for a by their large volume limits 

a = ((-d q V + 0- P ) P /(p 2 ) P (2') 



or 



a 



(t-p) P /(p 2 ) P (4') 



In general p depends on time, and so does a given by (2') or (4'). Suppose now that p is 
invariant under the IK or IE time evolution; then a and also V are time independent, so 
that 

0=(V) p = (d q V-q) p = Ud q V-p) p 

fit 

But then (2') and (4') coincide: the infinite system IK and IE evolutions have the same 
time invariant space ergodic states p. (Apart from the use of space ergodicity for an 
actually infinite system, this is the remark of Cohen and Rondoni [2]). 

Note that, if we replace in (3) m by m and V by V, imposing H = yields 

a =(-d q V-d q V + 0-^ 
m p z 

and in the infinite system limit we have again equivalence with the isokinetic ensemble. 
On the other hand, if H is not of the form p 2 /2m + V(q), the Gaussian thermostat doesn't 
give a term of the form — ap in (1) and we do not have equivalence with the isokinetic 
ensemble in the infinite system limit. 

For the purposes of nonequilibrium statistical mechanics one should presumably re- 
strict p to be an infinite system SRB state (defined so that the time entropy per unit 
volume is equal to the sum of the positive Lyapunov exponents per unit volume). Hope- 
fully, the space ergodic SRB states form a 2-parameter family parametrized by the average 
number of particles and the energy (or the kinetic energy) per unit volume. But the del- 
icate question of identifying the natural nonequilibrium steady states is here bypassed by 
the remark that they are the same for the infinite system IK and IE evolutions. 

In equilibrium statistical mechanics the proof of equivalence of ensembles is somewhat 
subtle, and uses in particular the concavity properties of the entropy (see [6]). One might 
think that the corresponding problem in nonequilibrium statistical mechanics would be 
even more difficult, and the above findings about the equivalence of IK and IE appear 
thus surprisingly cheap. What we have shown is however only that the IK and IE evo- 
lutions coincide (formally) in the infinite system limit; the detailed study of the natural 
nonequilibrium states remains to be made. 
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3. The constant a case. 

It is of interest to consider the equations (1) with a = constant. For this situation 
one obtains the following result. 

Proposition. 

Consider the evolution equations 

p = -d q V + £ - ap 
q = p/m 

in TM , where M C R" x and we impose elastic reflection on the boundary of M . We 
assume that a, m are constants > 0, and that V, £ are bounded. Then 

" msu P.-~(^ + V)< 1 ^T+ m a^ (5) 
2 max£ 2 

limsup f _ >00 p < 5 h 2m(max F — min V) (6) 

Furthermore, if the bounded measure p is invariant under time evolution and, and $ is 
any continuous function we have 



I 



p(dpdq)Q(£).(t-p-ap l ) = (7) 



From the evolution equations we obtain 

i(^ + V) = ^-p + d q -q=^-(-d q V + Z-ap) + d q -^ = ^-(Z-ap) (8) 
at 2m m m mm 

Let now e > and suppose that 



f- I V > tt— — r • n,axl - ( (9) 

2m 2ma z 



then 

p 2 > max( 2 /a 2 + e 

or 

a\p\ > max |£] + e' 
with e' > and thus, in view of (8), 

^(f + V)<-(m\-a\p\ 2 )<- l -^e' 
at 2m m m 



Therefore, as long as (9) holds, we have 



dV2m ' ~ 

for some S > 0, proving (5). From (5) we obtain immediately (6). 
Let ^' = $ then, by the ergodic theorem, 

J p(dp dg)*(|l).(£ ■ p - ap 2 ) = T lim I £ dt $(f^)-(£ • P ~ ap 2 ) 

= lim Z [ T dt*'(?-)4(^ + V)= lim £ f T dtU(f + V) 
T^ooTJo y 2m' dV2m ' t^oo T J dt v 2m ; 

2 

= lim ^[*(^- + y)]°° = o 

because \& is bounded in view of (5). This proves (7). [] 
Acknowledgements. 

I am indebted to Eddie Cohen, Giovanni Gallavotti, and Oscar Lanford for enlighten- 
ing discussions concerning this note. 

References. [1] N.I. Chernov, G.L. Eyink, J.L. Lebowitz, and Ya.G. Sinai. "Deriva- 
tion of Ohms law in a deterministic mechanical model." Phys. Rev. Letters 70, 2209- 
2212(1993). 

[2] E.G.D. Cohen and L. Rondoni. "Note on phase space contraction and entropy 
production in thermostatted Hamiltonian systems." Chaos 8,357-365(1998). 

[3] D.J. Evans and CP. Morriss. Statistical mechanics of nonequilibrium fluids. Aca- 
demic Press, New York, 1990. 

[4] G. Gallavotti. "Dynamical ensembles equivalence in fluid mechanics." Physica D 
105,163-184(1997). 

[5] G. Gallavotti. "Chaotic dynamics, fluctuations, non-equilibrium ensembles." 
Chaos 8,384-392(1998). 

[6] D. Ruelle. "Correlation functionals." J. Math. Phys. 6,201-220(1965). 

[7] Ya.G. Sinai. "A remark concerning the thermodynamic limit of the Lyapunov 
spectrum." Int. J. Bifurc. Chaos 6,1137-1142(1996). 



6 



